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Abstract 

We consider reshaping an obstacle virtually by using transformation optics 
in acoustic and electromagnetic scattering. Among the general virtual reshaping 
results, the virtual minification and virtual magnification are particularly studied. 
Stability estimates are derived for scattering amplitude in terms of the diameter of 
a small obstacle, which implies that the limiting case for minification corresponds 
to a perfect cloaking, i.e., the obstacle is invisible to detection. 

1 Introduction 

Since the pioneering work on transformation optics and cloaking [ffi \10\ [T5| 124], there 
is an avalanche of study on designs of various striking cloaking devices; e.g., invis- 
ibility cloaking devices [HJ Q2] ; field rotators [T] ; concentrators [T7] ; electromagnetic 
wormholes [H(5]; superscatterers [25], etc.. We refer to a most recent survey paper 
[7] for a comprehensive review and related literature. The crucial observation is that 
certain PDEs governing the wave phenomena are form-invariant under transforma- 
tions, e.g., Hemholtz equation for acoustic scattering and Maxwell's equations for 
electromagnetic scattering. Hence, one could form new acoustic or EM material pa- 
rameters (in the physical space) by pushing forward old ones (in the virtual space) via 
a mapping F. Such materials/media are called transformation media |24j . It turns 
out that the wave solutions in the virtual space with the old material parameters 
and in the physical space with the new material parameters are also related by the 
push-forward F. Those key ingredients pave the way for the design of optical devices 
with customized effects on wave propagation. 

In this paper, we shall be concerned with cloaking devices for acoustic and electro- 
magnetic obstacle scattering. As is known, there are two types of scatterers which are 
under wide study for acoustic and electromagnetic scattering, namely, the penetrable 
medium and the impenetrable obstacle. For a medium, the acoustic or EM wave can 
penetrate inside, and basically the medium accounts for the coefficients in the govern- 
ing PDEs. Whereas for an obstacle, the acoustic or EM wave cannot penetrate inside 
and only exists in the exterior of the object, and the obstacle is related to the domain 
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of definitions for the governing PDEs. The cloakings for acoustic or EM media have 
been extensively studied in transformation optics in existing literature and the theory 
has been well-established, we again refer to the review paper [7] for related discussion. 
For our current study, the cloakings for obstacles are considered and it is shown that 
the domain of definitions for certain PDEs can also be pushed forward under trans- 
formations. Using the transformation optics, one can push forward an obstacle in the 
virtual space to form a different obstacle in the physical space, and the ambient space 
around the virtual obstacle is then pushed forward to a cloaking medium around the 
physical obstacle. With a suitable push-forward F, it is shown that the scattering 
amplitude in the physical space coincides with that in the virtual space. That is, if 
one intends to recover the physical obstacle after being cloaked by the corresponding 
scattering measurements, then the reconstruction will give the image of the obstacle 
in the virtual space, but not the physical one, namely, the physical obstacle is virtu- 
ally reshaped with the cloaking. Principally, it has been shown that one can achieve 
any desired virtual reshaping effect provided an appropriate transformation F can be 
found between the virtual space and the physical space. 

Particularly, we consider virtually magnifying and minifying an obstacle. By 
magnification, we mean that the size of the virtual obstacle is larger than that of the 
underlying physical one. That is, under acoustic and EM wave detection, the cloaking 
makes the obstacle look bigger than its original size. Whereas by nullification, we 
actually mean virtually shrinking the obstacle, that is, the size of the virtual obstacle 
is smaller than that of the physical one. In the limiting case of minification, the 
virtual obstacle collapses to a single point, and this formally corresponds to a perfect 
cloaking, namely, the physical obstacle becomes invisible to detection. We note that 
in this case, the push-forward F blows up a single point in the virtual space to a 
'hole' (which actually is the physical obstacle) in the physical space. Hence, the map 
F is intrinsically singular, and the obtained transformation medium is inevitably 
singular. Correspondingly, the transformed PDEs in the physical space are no longer 
uniformly elliptic which also becomes singular. Therefore, in order to rigorously justify 
the perfect cloaking, we need to deal with the singular PDEs. Basically, one would 
encounter the same problems in treating perfect cloakings for acoustic or EM medium 
and several approaches are proposed to deal with such singularities. For perfect 
cloaking of conductivity equation, which can be considered as optics at zero frequency, 
the invisibility is mathematically justified in [10] by using the removability of point 
singularities for harmonic functions; whereas an alternative treating is provided in 
[12j . where near-invisibility is introduced from a regularization viewpoint and the 
invisibility is rigorously justified based on certain stability estimates for conductivity 
equation with small inclusions. For the finite frequency novel notion of 

finite energy solutions is introduced in [8] and the invisibility cloaking of acoustic 
and electromagnetic medium are then justified directly. For the perfect cloaking 
of obstacles considered in the present paper, we shall follow the approach in [12] to 
mathematically justify the invisibility by taking limit of near-invisibility. To that end, 
we derive certain stability estimates for scattering amplitudes in terms of the diameter 
of a small obstacle in both acoustic and EM scattering. Those stability estimates are 



2 



then used to show that the limiting process of minification cloaking corresponds to 
a process of near-invisibility cloaking, which in turn implies the desired invisibility 
result of the perfect cloaking. For practical considerations, all our reshaping studies 
are conducted within multiple scattering, that is, there is more than one obstacle 
component included. 

Finally, we would like to mention some unique determination results in inverse 
obstacle scattering, where one utilizes acoustic or electromagnetic scattering measure- 
ments to identify an unknown/inaccessible obstacle. The uniqueness/identifiability 
results correspond to circumstances under which one cannot virtually reshape an ob- 
stacle. In the case that the obstacle is situated in a homogeneous background medium, 
the uniqueness theory for inverse obstacle scattering is relatively well established, 
and we refer to |16j for a survey and relevant literature. Whereas in |llj.[15J.|22j. 
the recovery of an obstacle included in certain inhomogeneous (isotropic) medium is 
considered. It is shown in [11] and [22] that if the isotropic medium is known a priori, 
then the included obstacle is uniquely determined by the associated scattering am- 
plitude. Under the assumption that the isotropic medium and the included obstacle 
has only planar contacts, it is proved in [15] that one can recover both the medium 
and the obstacle by the associated scattering amplitude. The argument in [15] also 
implies that an obstacle surrounded by an isotropic medium cannot produce the same 
scattering amplitude as another pure obstacle. This result essentially indicate that 
transformation media for virtually reshaping an obstacle must be anisotropic. 

The rest of the paper is organized as follows. In Section we consider the re- 
shaping for acoustic scattering, where virtual minification and magnification are first 
considered consecutively, and then we present a general reshaping result. Similar 
study has been conducted for reshaping a EM perfectly conducting obstacle in Sec- 
tion on 



2 Virtual reshaping for acoustic scattering 
2.1 The Helmholtz equation 

Let M be an open subset of M 3 with Lipschitz continuous boundary dM and con- 
nected complement M + := ]R 3 \M. Let (M + ,g) be a Riemannian manifold such that 
g is Euclidean outside of a sufficiently large ball Br containing M. Here and in the 
following, Br shall denote an Euclidean ball centered at origin and of radius R. In 
wave scattering, M denotes an impenetrable obstacle and the Remannian metric g 
corresponds to the surrounding medium with the Euclidean metric go := 5j repre- 
senting the vacuum. In acoustic scattering, a = (c 4 -') 3 , )=1 with er ij := \f\g\g 1 -' is the 
anisotropic acoustic density and y/\g\ = \<r\ is the bulk modulus, where (<7 l3 ")f j =1 is 
the matrix inverse of the matrix (<7y)f ,- = i, and \g\ = detg, \a\ = deter. Formally, we 
have the following one-to-one correspondence between a material parameter tensor 
and a Riemannian metric 

<7« = \g\ 1/2 g ij or g» = H~V^'. (2.1) 
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We consider the scattering for a time-harmonic plane incident wave u l = exp{ikx ■ 
9}, 8 £ § 2 due to the obstacle M together with the surrounding medium (M + ,g). 
The total wave field is governed by the Helmholtz equation 

A g u + k 2 u = Q in M+, (2.2) 
u\ dM = 0, (2.3) 

where the Laplace-Beltrami operator associated with g is given in local coordinates 

by 



The homogeneous Dirichlet boundary condition (|2.3p means that the wave pressure 
vanishes on the boundary of the obstacle. M is usually referred to as a sound-soft 
obstacle. The scattered wave field is as usual assumed to satisfy the Sommerfeld 
radiation condition. Taking advantage of the one-to-one correspondence (|2.ip be- 
tween (positive definite) acoustic densities a and Riemannian metrics g, we proceed 
to mention a few facts about the form-invariance of the Hemholtz equation under 
transformations. For a smooth diffeomorphism F := f^i — ► 0,2, V = F(x), the metric 
g{x) transforms as a covariant symmetric 2-tensor, 



(2-4) 



Q ryr> I Q rj* TTl 

~ 9lj {y) := (F*g)ij(y) = ^ g lm 

l,m=l U 9 

and then, for u = u o F, we have 

(A g + k 2 )u = (Ag + k 2 )u = 0. (2.5) 

Alternatively, using (|2.ip , one could work with the Helmholtz equation of the following 
form 

3 

di(a ij d jU ) + k 2 \a\u = 0, (2.6) 

and then, for u = u o F, we have 

3 

di(a ij d jU ) + k 2 \a\u = 0. (2.7) 

Here cr is the push- forward of a which, by using (|2.1|) and (|2.4|) . can be readily shown 
to be given by 

where DF denotes the (matrix) differential of F and (DF) T its transpose. 

Throughout, we shall work with a £ L°°(M + ) 3x3 and F is orientation-preserving, 
invertible with both F and i 7-1 (uniformly) Lipschitz continuous over M + . So, it is 
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appropriate to work with the following Sobolev space for the scattering solution to 

d23D-dl3D, 



H} oc (M + ) = {u£ @'(M + );u G H l (M + n B p ) for each finite p with M C B p }. 

The system (j2.2[) - ()2.3|) . or (|2.6p and (|2,3p is well-posed and has a unique solution 
u G Hl oc (M + ) (see [H]). Noting that the corresponding metric outside a ball -Br D M 
is Euclidean, we know u is smooth outside Br. Furthermore, the solution u(x, k, 9) 
admits asymptotically as |x| — ► +00 the development (see [3]) 

u(x, 9, k) = e ikxS + ^-A(0', 9, k) + O(^), (2.9) 
\x\ \x\ z 

where 0' = x/\x\ G S 2 . The analytic function A(9',6,k) is known as the scattering 
amplitude or far-field pattern. According to the celebrated Rellich's theorem, there 
is a one-to-one correspondence between the scattering amplitude A(9', 9, k) and the 
wave solution u(x,6,k). Throughout, we consider the scattering amplitude for the 
virtual reshaping effects. 

We shall denote by M Q)(M + , g) a cloaking device with an obstacle M and the 
corresponding cloaking medium (M + ,g). The metric g is always assumed to be 
Euclidean outside a sufficiently large ball containing M , namely, the cloaking medium 
is compactly supported. If we know the support of the cloaking medium, say M'\M, 
we also write M Q)(M'\M, g) to denote the cloaking device. 

Definition 2.1. We say that M Q)(M + , g) (virtually) reshapes the obstacle M to 
another obstacle M, if the scattering amplitudes coincide for M ®(M + , g) and M, 
i.e. 

A{9\ 9, k; M 0(M+, 5 )) = A(9', 9, k; M). 

We would like to remark that according to the correspondence (]2.ip . the cloak- 
ing device in Definition (|2.ip can also be written as M Q)(M + ,o~), where a is the 
(anisotropic) acoustic density for the cloaking medium. 

2.2 Virtual minification by cloaking 

We first consider the reshaping effects for a special class of obstacles, which are star- 
shaped and referred to as l p -ball shaped obstacles in the following. They are domains 
in M 3 of the form 

{x G M 3 ; \\x\\ p = r}, 
where p G [1, +00], r > is a constant and for x = (xi,X2,x^) 

( 3 

ll^llp — ( ^ ^ \%i I 
\i=l 

Obviously, || • H2 = | • | and an Z 2 -ball is exactly an Euclidean ball. For Z p -ball shaped 
obstacles, we can give the transformation rule explicitly and correspondingly, the 
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Figure 1: Illustration for nullification: the cloaking device B Rl ©(-B^ , g) virtually 
reshapes B Rl to -B_r . 

cloaking material parameters for those obstacles can be derived explicitly. Henceforth, 
we write B R)P to denote an / p -ball of radius R and centered at origin, whereas as 
prescribed earlier, we write Br := Br 2- We also denote by M + the complement of a 
domain M in M 3 . 

Let M = B Rup with Ri > 0. Let Rq,R 2 be such that < R < R\ < R 2 . We 
define the map, F : B^ op i-> B~^ ip by 

_ „ f v _ / 2/, ^ \\y\\ p > R 2 , 

X ''~ {y) ~\(^ Rz + ^\\y\\ P )A> torR < Mp <R 2 . (2 - 10) 

It is noted that F is (uniformly) Lipschitz continuous over B^ Q and maps B R2iP \B RoiP 
to B R2tP \B Rup . For i?i < < R 2 , let 

gi (x) = (F*go)(x). (2.11) 

Set 

g(x) = \ 91{X) ' f ° rEl || <ll4<i?3 ' (2.12) 
|jo(4 for > i? 2 . 

We have 

Proposition 2.2. XTie cloaking device B RljP ^(B^ ,g) with g defined in 112.12)) . re- 
shapes B RltP virtually to B RlhP . That is, the physical obstacle B Rl ^ p with the cloaking 
material (B R2tP \B RltP , g\) is virtually minified to the obstacle B RotP with a minifica- 
tion ratio K := Rq/R\ (see Fig. U\for a schematic illustration). 

Proof. Let u{y) G H^J^B^ ) be the unique solution to the Hemholtz equation f)2 . 2|) - 
(|2.3p corresponding to the obstacle B Ro ^ p . Whereas, we let v{x) G H^^B^ ) be the 
unique solution to the Helmholtz equation (|2.2p - (|2.3p corresponding to B RljP ©(-B^i p , 
g). Define u(x),x € B^ , be such that u = u o F, i.e. u = F^u = (F~ 1 )*u. It is 
clear that u £ Hi oc (B^ ) since i 7 is bijective and both i* 1 and i 7-1 are (uniformly) 
Lipschitz continuous. Moreover, noting F(dB Ro , p ) = dB Rl)P , we know u\dB R p = 0. 

By the invariance of Helmholtz equation under transformation, it is readily seen 
that u = v. Hence, 

A(9',9,k-B RuP Q)(B+ up ,g))=A(6>,e,k;B Ro , p ). (2.13) 
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□ 

For an Euclidean ball B Ro C M 3 , by separation of variables, we have 



A(6', 9, k; B Ro ) = - V(2n + 1) J " 0) P ra (cosV), (2.14) 



where j n (t) and hn\t) are respectively, the n-th order spherical Bessel function and 
spherical Hankel function of first kind, P n (t) is the Legendre polynomial and ip = 
Z(9,9'). Using the asymptotical properties 

j n (t) = 0(t n ), h^(t) = 0(t- n - 1 ), n = 0,l,..., ast^+0, 

it is straightforward to show 

A(9',9,k;B Ro ) = 0(R ) as R - +0, (2.15) 

Now we consider the limiting case for minification, namely k — > +0 or equivalently 
#0 -> +0. By <|213|> and <|21K|> . 

9, k; B Rl 0(B+ , 9 )) = A(0', 0, fc; B^) = 0(i? o ) (2.16) 

as Rq — > +0. That is, 

Proposition 2.3. The limit for minification of an Euclidean ball B Rl in Proposi- 
tion \2.2\ gives a perfect cloaking, namely, it makes the obstacle invisible to detection. 

In the limiting case with k = 0, the transformation in f|2.10j) becomes 

{y, for ||y||p > R%, 
c^+^Wnit. fOT < \\v\\p < (2 " 17) 

which maps IR 3 \{0} to IR 3 \-Br 1iP , i.e., it blows up the single point {0} to B RljP . It is re- 
marked that the map H in (|2.17p with p = 2 is exactly the one used in [21 [TO] for perfect 
cloaking of conductivity equation, and in |24j for perfect cloaking of electromagnetic 
material tensors. Next, we take the case with p = 2 as an example for a simple analysis 
of the perfect cloaking medium. The corresponding metric (H*go)(x) in Br 2 \B Ri is 
singular near the cloaking interface, namely dB Rl . In fact, considering in the standard 
spherical coordinates on Br 2 \{0}, (r, </>,#) i— ► (r sin 9 cos (f>, r sin 9 sin cj>, r cos 9) £ M 3 
and by (12. 4ft . it can be easily calculated that 

/ A 2 
g : = H*go = A 2 (r - R x ) 2 

\ A 2 (r- 

where A = -R2A-R2 — Ri)- That is, g has one eigenvalue bounded from below (with 
eigenvector corresponding to the radial direction) and two eigenvalues of order (r — 
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R\) 2 approaching zero as r — ► +R\. Hence, if the perfect cloaking is analyzed directly, 
one needs to deal with the degenerated elliptic equation near the cloaking interface. 
So, a suitable choice of the class of weak solutions to the singular equation must be 
purposely introduced, as the finite energy solutions considered in [8] for invisibility 
cloaking devices of acoustic and electromagnetic media. Clearly, our earlier analysis 
on the perfect cloaking of an Euclidean ball avoid singular equation by taking limit. 
This is similar to [12] for the analysis of perfect cloaking of conductivities in electrical 
impedance tomography by regularization. Here we would like to point out that there 
is no theoretical result available showing that the limit of the regularized solutions 
obtained by the approach of the current paper by sending n — > are the finite energy 
solutions in the sense of [8] . A further study in this aspect may provide more insights 
into the invisibility cloaking. 

In order to achieve the similar invisibility result for a general / p -ball shaped ob- 
stacles, we need to derive stability estimates similar to (12.160 for generally shaped 
obstacles with small diameters. This is given by Lemma 12.41 below, proved using 
boundary integral representation rather than separation of variables, and the obsta- 
cles could be generally star-shaped. On the other hand, from a practical viewpoint, 
we consider the scattering with multiple scattering components and only some of the 
components are cloaked. We shall show that the virtual reshaping takes effect only 
for those cloaked components and the other uncloaked components remain unaffected. 
Particularly, those perfectly cloaked components will be invisible, even though there 
is scattering interaction between the obstacle components. We are now in a position 
to present the key lemma. In the sequel, we let B be a simply connected set in R 3 
whose boundary is star-shaped with respect to the origin of the form <9B = 5ro(9)9, 
where 9 G S 2 ,r (6») G C 2 (§ 2 ) and 6 > 0. Let B be the domain {x G K 3 ; \x\ < r (9)}. 

Lemma 2.4. Suppose M\ n Bo = 0, then we have 



Proof. Let y) = e ' x y ' / (4tt\x— y\) be the fundamental solution to the Helmholtz 



operator (A + k 2 ). We know that u(x; k, 9) G C 2 (M 3 \Afi U B) n C(lR 3 \(Mi U B)) and 
can be represented in the form (see [3]) 



where ipi ^ C(dM\) and <f2 G C(<9B) are density functions, and rj ^ is a real 
coupling parameter. The densities (p± and if2 are unique solutions to the following 
integral equation (see [3]) 



A(9',9,k;M 1 iJM) = A(6 , ,e,k;M 1 ) + 0(6) as 5 -> +0. 



(2.18) 




(2.19) 




(2.20) 



+2 




dv{y) 



irj<&{x,y) >ip 2 (y) ds(y) 



= -2e 



ikx-8 
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for x G dM\ U &B, where ip(x) := yi(^) for x G <9Mi and y(a;) := ^2^) fo r £ €E <9B. 
We introduce the integral operators 

(S m )(x) = 2 / ^(x,y)My)ds(y), (K Wl )(x) = 2 / ^2A^{y)ds{y) 

JdMi JdMi av \V) 

(S 2 <p 2 )(x) = 2 [ *(x,y)<p2(v)Mv)> (K 2 cp 2 )(x) = 2 [ ^^ 2 (y)<fe(y), 
Jan Job 9u{y) 

and set 

h^x) := -2e ikx - e , x G dM x \ h 2 (x) := -2e ikx -\ x G 01. 

Then equation (|2.20p can be rewritten as 

[<pi + Kiipi - iS\ipi + K 2 f 2 - i7]S 2 ip 2 ]{x) = hi(x), x G dMi ( 2 -21) 
[ip 2 + K 2 ( P2 -ir]S 2 (p 2 + K 1 (p 1 -iS 1 (p 1 ](x) = h 2 (x), x G dM. (2.22) 

It is remarked that the integral operators involved in equations (12.21j) and (|2.22l) 
with weakly singular integral kernels have to be understood in the sense of Cauchy 
principle values and we refer to [3] and [19] for related mapping properties. Clearly, 
(pi and ip 2 are functions dependent on 5. We next study their asymptotic behaviors 
as 5 — > +0. To this end, we fix 5 > but being sufficiently small and take rj = <5 _1 . 

In the sequel, without loss of generality, we may assume that dist (dM±, c?Bq) > 
Co > 0, otherwise one can shrink Bo to l/2Bo- By straightforward calculations, it can 
be easily shown that 

\\K 2 - ^52||c(aB)-^C(3Afi) = 11-^1 ~~ i Si\\ C (dM 1 )-*C(dM) = 0(1). (2.23) 

Next, for x G <9Bo, we define 

(K <f>)(x) = 2 [ d ^\ V h (y) ds(y), (S <P)(x) = 2 [ $ (x, y)<f>(y) ds(y), 
JdMo dv(y) J dBo 

where <p 6 C(8Mq) and &o(x,y) = 1/(4tt\x — y\) is the fundamental solution to the 
Laplace operator. It is known that both So and Kq are compact operators in C(cffio) 
(see [2]). By changing the integration to the boundary of the reference obstacle <9Bo, 
we have 

p ik\x—y\ r ik5\x'-y'\ 

(S 2 <p)(x) =2 / My) ds(y) = 25 / — -<p(6y') ds(y'), 

Jan F - v\ JdMo F - y \ 

(k 2 <p)(x) =2 / a[ , /MvMv) My) 

Jan \ F - y\ J 

( ikS\x'-y'\\ 
d -n jr /dv(y>M5y>) ds(y'l 

, V f - y \ ) 
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where ip G C(dM), x,y G <9B and x' := x/5,y' = y/5 G 8Mq. Then by using power 
series expansion of the exponential function e ikS \ x ~y\ } we have by direct calculations 

||-52 - S \\c(8dM )~*C{dM ) = C((5), \\K 2 - K \\c(SdM )^C(dB ) = 0(5 2 ). (2.24) 

By changing the integration to <9!>o and using the results in (|2.24p . we have from 
(12T22D that 

ip 2 (5x) = [I + K Q -iS + O(5)}- 1 [h 2 (5x)-(K 1 -iS 1 )ip 1 (5x)], x G 8M Q . (2.25) 

It is noted here that (I + Kq — iSo) is bounded invertible (see [2]). Then, plugging 
(12.25ft into (I2.21|) and using the relations in (|2,23ft . we further have 

<pi = [I + K x - iSi + OiS)]- 1 ^! + 0(S)), (2.26) 

which, by noting I + K\ — \S\ is invertible (see [3]), gives 

<pi = <pi + 0(S), (2.27) 

where 

ft = [I + Kx-iS^hx. 
Furthermore, ([2.27ft together with ([2.25ft implies that 

W = 0(1). (2.28) 

Finally, by (|2.19p we know 

1 f f de- iky - 9 ' 



iff f)p~ y 1 



(2.29) 



Using the estimates in (j2.27ft and ([2.28ft to (|2.29ft and changing the integration over 
<9B to <91B>o, we have 

A(8', 9, k; Mi UB) = A(0', 9, k; Mi) + 0(5), 
where we have made use of the fact that 

A(9>,9,k;M 1 ) = - / ^(y)\^—^- te -^'y\ds{y). 

4tt J dMl { du(y) J 

The proof is completed. □ 

Remark 2.5. If dM\UdM is only Lipschitz continuous (whence r$(0) G C 0,1 (S 2 ) ), one 
can make use of the mappiny properties of relevant boundary layer potential operators 
presented in /TP) / and derive similar estimate. 
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Br B Hl 

Figure 2: Illustration for magnification: the cloaking device B Ro ®(-Bro,<?) virtually 
reshapes B Ro to B Rl . 

Proposition 2.6. Suppose that M\ n B R2tP = ^ e cloaking device ((Mi U B Rl ^ p ) 
0((Ml U B RuP ) + ,g)), where g(x) is gi in V2.11\) for R 1 < \\x\\ p < R 2 and g (x) for 
x £ (Mi U Bji 2iP ) + , reshapes the obstacle Mi U B RltP to M±U Bj^^. Furthermore, the 
limiting case with Rq = corresponds to the perfect cloaking of B Rl)P , namely 

A(9', 9, k- (Mi U B RljP ) 0((Mi U B RuP ) + ,g)) = A(9', 9, k; Mi). (2.30) 

Proof. Let F be the transformation in (|2.10p and let F be the restriction of F over 
(M\\JB Ro ) + . Clearly, g = F*go. By a similar argument as the proof of Proposition ^. 21 
it is easily seen that Mi U B RltP is virtually reshaped to Mi U B Ro;P by the cloaking 
of g, i.e., 

A(9', 9, k; (Mi U B RltP ) 0((M a U B RljP ) + ,g)) = A(6', 9, k; M l U B Ro , p ). 

Next, by Lemma 12.41 

A(9',9,k;M 1 \J B^p) = A(e',6,k;Mi) + 0(Ro) as R +0, (2.31) 

and hence the limiting case with Rq = yields an ideal cloaking of B RltP . □ 

Remark 2.7. Clearly, Proposition \2. 61 implies a same invisibility result for perfectly 
cloaking an l p -ball as that in Proposition^^ for perfectly cloaking an Euclidean ball. 

2.3 Virtual magnification by cloaking 

Let < Ro < Ri < i?2, and let B R(I)P be the obstacle which we intend to virtually 
magnify to B RltP by using a cloaking for B Ro ^ p supported in B^^Bj^^ (see Fig. [2]). 
We define r = Ri/Ro to be the magnification ratio. 
Let K : R 3 \B RuP ■-> R 3 ^^ be defined by 



x := 



V, for ||y|| p > i? 2 , 

& 1 ** + %F%Wv\\p)A> for A)< ||y|| p < i* 2 . 



n 



It is verified directly that K maps the P-annulus R\ < \\y\\ p < R2 to the P-annulus 
Ro < \\y\\p R2- Moreover, K is bijective and both K and K -1 are (uniformly) 
Lipschitz continuous. Set 

I = K*g , (2.32) 
be the metric in 2?^ . Clearly, g is Euclidean outside Br 2)P . 

Proposition 2.8. The cloaking device Bji (hP Q)(Bii 2tP \BR (hP , g) with g defined in 
h2.'3ty) . reshapes Br 0iP virtually to Br 1)P . That is, the physical obstacle Bj^ tP with the 
cloaking material (Bn 2tP \Bji (hP , g) is virtually magnified to the obstacle -Br 1iP with a 
magnification ratio r := R\/Rq. 

Proof. Let u G H} oc ( b r 1)P ) be tne solution to the Helmholtz equation §F2$-§FS\) 
associated with the obstacle Br 1)P . Define u = K*u € b I O c(Rr 0P )- Again, by the 
invariance of the Helmholtz equation under transformation together with the fact 
that u\ob r p = u\qb r = 0, we see u is the scattering solution corresponds to 

3*0* ©(^i). That is, 

A(9', 9, k; B Ro , p 0(B+ p , $)) = A(9', 9, k; B RuP ). 

□ 

In Proposition ^. 2\ we use F to compress the vacuum to achieve a transformation- 
based minification device, whereas in Proposition 12.81 we use K to loosen up the 
vacuum to achieve a transformation-based magnification device. Note that R2 > 
R±, the cloaking device is of size larger than the virtual obstacle image, though 
the virtual obstacle could be of size arbitrarily close to the cloaking device. Hence, 
the cloaking in Proposition 12.81 is not of magnification in the real sense. However, 
our magnification result is still of particular practical interests, e.g., if one is only 
interested in recovering an obstacle without knowing a priori that it is cloaked, then 
the scattering reconstruction will give a virtually magnified obstacle. On the other 
hand, we would like to mention that in [20l EU [23] , it is demonstrated that a coated 
cylindrical core can be extended beyond the cloaking shell into the matrix, where 
the cloaking material must be negative refractive indexed, namely, the corresponding 
metric g has negative eigenvalues. A general strategy is presented in fl4J on how 
to devise a negative refractive indexed (NRI) cloaking by using the transformation 
optics. There, the transformation F is neither injective nor orientation-preserving, 
which maps a right-handed medium to left-handed medium. Based on NRI cloaking, 
it is shown in |18} [25] that one can virtually reshape a cylindrical perfect conductor 
of size bigger than the cloaking device. However, all the aforementioned results are 
essentially based on exerting transformation directly to the analytical solutions, which 
is not of the main theme of the present paper. 

2.4 Virtually reshaping acoustic obstacles by cloaking 

Our discussion so far has been mainly concerned with the minification and magnifica- 
tion of obstacles by cloaking. Clearly, we may consider virtually reshaping an obstacle 
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arbitrarily provided a suitable transform can be found with which we can make es- 
sential use of the transformation invariance of the Helmholtz equation. Let M be an 
obstacle with m pairwise disjoint simply connected components Mi, I = 1,2, ... , m, 
i.e., M = M h Let M z C M{ and M{ n M' v = 0, for 1,1' = 1,2, ... ,m and I / Z'. 
Set M' = U™i M/. Let M = (Ji=i M be another obstacle with M L C Afy. Suppose 
there exist 

Fi : M 3 \M, i — y M 3 \M Z , Z = 1, 2, . . . , m 

such that Fi is orientation-preserving and invertible with Fi and -F^ 1 Lipschitz con- 
tinuous, and Fi = id outside M[. Set g\ = (Fi)*go and let Mj 0(M/\M/, <#) be a 
cloaking device for Afj. 

Theorem 2.9. TTie cloaking device M' virtually reshapes the obstacle M to M . That 
is, 

A{9', 6, k; M 0(M'\M, g')) = A(9', 6, k; M), 

where g' is g\ in M[\M\. 

The proof is already clear from our earlier discussion on minification and magni- 
fication. We have several important consequences of the theorem. 

Remark 2.10. Suppose that some of the components of M are uncloaked, say Mi 
for 1 < I < m' < m, and this corresponds to taking Mi = M[ = M\ and Fi = id for 
I = 1,2,.. . ,m'. 

Remark 2.11. If for some Mi being star-shaped w.r.t. certain point, and the trans- 
formation Ff 1 shrinks Mi only in the radial direction to Mi, then the case with M 
degenerated to a single point corresponds to an ideal cloaking for Mi . By using a sim- 
ilar argument as that for Proposition \2.6\ together with the estimate in Lemma \2.4\ 
one has that Mi is invisible to detection. In fact, by repeating the argument, the 
same conclusion holds when there are more than one obstacle component is perfectly 
cloaked. 

It is noted that in Remark I2.11|, the perfectly cloaked obstacle components are 
required to be star-shaped, and this is because we need to make use of the estimate 
in Lemma 12.41 to achieve the invisibility. In order to show the prefect cloakings of 
more generally shaped obstacles, one may need different thoughts. 

In the rest of this section, we shall indicate that all our previous results on vir- 
tual reshaping in space dimension three can be straightforwardly extended to the 
two dimensional case. In fact, for two dimensional scattering problem, the (positive 
definite) acoustic density a £ L°°(Af + ) 2x2 also transforms according to (|2.8p . There- 
fore, the reshaping result presented in Theorem 12.91 is still valid in M 2 . In order to 
achieve invisibility for perfect cloaking of star-shaped obstacles in M 2 , one needs to 
show a similar estimate to Lemma 12.41 Indeed, replacing by the first kind 

Hankel function jH^\k\x — y\) of order zero in the proof of Lemma 12.41 and using 
the corresponding mapping properties of the integral operators involved (see [3]), one 
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can obtain by similar arguments the following estimate to the scattering problem in 
M 2 (see (|2.18p for comparison), 

A(9',e,k;Mx UB) = A(0' , 6 , k; Mx) + 0{\log5\- 1 ) as S -» +0. (2.33) 

Obviously, with (|2.33[) one can show that the perfect cloaking of a star-shaped obsta- 
cles in M 2 makes it invisible to detection. 



3 Virtual reshaping for electromagnetic scattering 
3.1 The Maxwell's equations 

We define Maxwell's equations for the scatterer M @(M + ,g) as the one introduced 
in Section \2. 11 Using the metric g, we define a (positive definite) electric permittivity 
tensor e and magnetic permeability tensor il by 

e*i = LL ij = \g\ 1/2 g ij on M+. (3.1) 

It is clear that e = j=i an d A* = (/%)fj=i are invariantly defined and transform 
as a product of a (+)-density and a contravariant symmetric two-tensor with the 
same rule as that for acoustic density a in (|2.8p . We consider the scattering due to 
the scatterer M Q(M + , g) corresponding to some incident wave field. The resulting 
total electric and magnetic fields, E and H in M + , are defined as differential 1-forms, 
given in some local coordinates by 

E = Ej dx j , H = Hj dx j . 

Here and in the following, we use Einstein's summation convention, summing over 
indices appearing both as sub- and super-indices in formulae. Then (E, H) satisfies 
Maxwell's equations on (M + , g) at frequency k 

dE = ik * g H, dH = -ik * g E, (3.2) 

where * g denote the Hodge-operator on 1-forms given by 

* g (Ej dx j ) = ^\g\ 1/2 g jl EjSi pq dx p A dx q = ^£ jl EjSi pq dx p A dx q , 

with sipg denoting the Levi-Civita permutation symbol, and si pq = 1 (resp. si pq = — 1) 
if (l,p,q) is an even (resp. odd) permutation of (1,2,3) and zero otherwise. By 
introducing, for H = Hjdx J , the notation 



the exterior derivative may then be written as 

f) J-f 1 

d{H q dx q ) = -g^dx p A dx q = -(curl H) l si pq dx p A dx q . 
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Hence, in a fix coordinate, the Maxwell's equations f|3.2|) can be written as 

(curl E) 1 = ikn jl Hj, (curl H) 1 = -ike jl Ej. (3.3) 

Without loss of generality, we take the incident fields to be the normalized time- 
harmonic electromagnetic plane waves, 

EHx) := yCVLvl cmlpe lkxe , H'Hx) : = cm\pe ikxe , 
k 

where p G M 3 is a polarization. As usual, the radiation fields are assumed to satisfy 
the Silver-Miiller radiation condition. To complete the description, we further assume 
that the obstacle M is perfectly conducting, and we have the following two types of 
boundary conditions on dM: the perfect electric conductor (PEC) boundary condition 

v x E\ 9M = 0, 

or the perfect magnetic conductor (PMC) boundary condition 

v x H\ 9M = 0, 

where v is the Euclidean normal vector of dM. 

We shall work with e, \i G (L oc (M + )) 3x3 . It is convenient to introduce the follow- 
ing Sobolev spaces 

£T(curl; !!) = {«£ £ 2 (^) 3 ; curlu G L 2 (0)}, 
flioc(curl; M + ) = {u G ^'(M+) 3 ; u G F(curl; M + n B p ) 

for each finite p with M C B p }. 

Then it is known that there exists a unique solution {E,H) G # ioc (curl; M+) © 
7/; oc (curl; M + ) to the electromagnetic scattering problem. Moreover, the solution 
E(x,k,p,6) admits asymptotically as \x\ — > +oo the development (see [3]) 

E(x,k,p,9) =E i (x) + ^E 00 (e',k, P ,9) + 0( T1 =), (3.4) 

where 9' = x/\x\ G S 2 . The analytic function E^iO 1 \k,p,0) is known as the electric 
far-field pattern. Similar to Definition 12.14 we introduce 

Definition 3.1. We say that M Q)(M + , g) (virtually) reshapes the obstacle M to 
another obstacle M, if the electric far-field patterns coincide for M ®(M + , g) and 
M, i.e. 

E^e', k,p, 9- M 0(M+ g)) = Eootf, k,p, 9- M). 

We would also like to remark that the cloaking device in Definition 13.11 can also 
be written as M0(M + ,e,/i) according to the correspondence (|3.ip . where e and 
jj, are respectively, electric permittivity and magnetic permeability for the cloaking 
medium. 
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3.2 Virtually reshaping electromagnetic obstacles by cloaking 

We consider the virtual reshaping for electromagnetic obstacles by cloaking. Let M = 
UlLi Mi,M' = (JfcLi M[,M = U|2=i Mi and F u I = 1, 2, . . . , m be those introduced in 
Section 12.41 Furthermore, we assume that i 7 ) is normal-preserving in the sense that 

v l = v l oF h l = l,2,...,m, 

where ^ and are, respectively, the Euclidean normals to dM\ and <9M;. E.g., if M\ 
and M/ are both star-shaped w.r.t. the origin, say dM\ = f{9)9 and dM[ = r(9)9 
with f/r = c being some constant, then Fi is normal-preserving since one has 

r{6)9 - Gradr 

VI If [9)9 ~ Vl\r(9)9 - —=====. 

y r z + |Gradr| z 

Particularly, if is / p -ball shaped, the transformation of the following form 

F(y) = (a + b\\y\\ p )/ r , 

\\y\\p 

is normal-preserving, which transforms an P-ball of radius f into another P-ball of 
radius r = a + bf. 

Concerning the virtual reshaping, we have 

Theorem 3.2. The cloaking device M' virtually reshapes the obstacle M to M . That 
is, 

E^O', k,p, 9; M 0(M'\M, g 1 )) = E^O', k,p, 9; M), 
where g' is gi in M/\M;. 

Proof. Let F : M 3 \M i-> R 3 \M be such that F\ M ,^ = Fi\ M ,^,l = 1,2,... ,m 

and F = id over M 3 \M'. Let {E,H) G flj oc (curl; M+) © i?j oc (curl; M+) be the 
unique scattering solution corresponding to the perfect conducting obstacle M. Define 
E = F^E and H = F*H. Clearly, (E, H) £ F Zoc (curl; M+)©F /oc (curl; M+) according 
to our requirements on the mappings Fi's, I = 1,2, ... ,m. Moreover, noting Fi's, 
1 < I < m, are normal-preserving, we know v x E\qm = (resp. v x H\qm = 0) 
if M is a perfectly electric conducting obstacle (resp. perfectly magnetic conducting 
obstacle). Hence, (E,H) is the unique solution corresponding to the cloaking device 
M0(M'\M,g')- Therefore, we have 

E 00 (6', k,p, 9; M (M'\M, g')) = E oo (0', k,p, 9; M). 

□ 

With Theorem [321 all the virtual minification and magnification results for acous- 
tic obstacle scattering can be straightforwardly extended to the electromagnetic obsta- 
cle scattering. In order to obtain similar invisibility results for a perfectly conducting 
obstacle when some of its star-shaped components are perfectly cloaked, we need a 
lemma similar to Lemma 12.41 in the following for electromagnetic scattering. 
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Lemma 3.3. Let M\, Bo and B be the same as those in Lemma \2.4\ then we have 

£oc(#', k,p,0; Mi UB) = E oo (0',k,p,e;Mi) + 0{5) as 5 -> +0. (3.5) 

Proof. We first introduce the space T°' a (dMi),0 < a < 1, consisting of the uni- 
formly Holder continuous tangential fields a equipped with the Holder norm, and 
T% a {dMx) = {a E T°> a (dMi);DWa E C°' a (<9Afi)}. Similarly, one can introduce 
T°' Q (5Mi U<9B). We know the solution (£, if) E C°> Q (M 3 \(Afi UB)) ©C°> a (R 3 \(Mi U 
and can be expressed (cf. [3]) 

E(x) = E\x) + curl / ai(y)*(ac, y)ds(y) + i / x (Sga)(y)*(x, y)ds(y) 

+ curl / a 2 (y)$(x,y)ds(y) + in / i/(y) x (S$a 2 ){y)<Z>(x,y)ds(y) 



and if (a;) = curlf£(a;)/iA;, where rj ^ is a real coupling parameter. Here, So is 
the operator as defined in the proof of Lemma 12.41 but with the integration domains 
changed according to the context and the densities a\ E T®' a (dMi),a 2 E T® ,a (dM) 
satisfy 

(a + Mxax + iNiPSlai + M 2 a 2 + ir}N 2 PS%a 2 )(x) = V(x), x E dM\ U <9B (3.6) 

where a{x) := a\{x) for x E dM\ and a{x) := <22(x) for x E 9B, and V{x) := 
-2v x for PEC obstacle and V{x) := -2v x H\x) for PMC obstacle. The 

operators involved in (13. 6ft are respectively given by 



(Miai)(x) :=2 / z/(z) x curl x {ai(y)$(x, y)} ds(y), 

(A/i6i)(x) :=2i/(x) x curl curl / u(y) x 6i(y)$(x,y) ds(y), 
(-M 2 a 2 )(x) :=2 / x cml x {a 2 (y)$(x,y)} ds(y), 



(A/262) (x) :=2i/(x) x curl curl / v(y) x b 2 (y)$(x,y) ds(y), 

JdM 

Pc : = (v x c) x c. 

We again refer to 0[3] for relevant mapping properties of the above operators. Finally, 
a similar asymptotic analysis to that implemented in the proof of Lemma 12. 4[ one 
can complete the proof. □ 

Clearly, with Lemma 13.31 we have similar invisibility result for electromagnetic 
scattering as those remarked in Remark 12. Ill for acoustic scattering. 
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